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Effects  of  mechanical  compression  and  PTFE  content  on  the  through-plane  gas  permeability  of  gas  dif¬ 
fusion  layers  (GDLs)  of  PEM  fuel  cells  are  investigated  both  experimentally  and  theoretically.  A  new 
test  bed  is  designed  and  built  which  allows  pressure  drop  and  air  flow  rate  measurement  for  various 
GDL  samples.  The  measured  values  are  used  to  calculate  the  through-plane  permeability.  Various  GDLs 
are  obtained  and  tested  over  a  wide  range  of  PTFE  content  and  compression  ratio.  The  experimental 
data  show  a  reverse  relationship  between  the  through-plane  permeability  and  both  PTFE  content  and 
mechanical  compression.  An  existing  model  for  through-plane  permeability  of  planar  fibrous  structures 
is  revisited  to  develop  a  model  that  accommodates  effects  of  PTFE  content  and  mechanical  compression. 
The  proposed  model  captures  the  trends  of  the  experimental  data  for  the  through-plane  permeability, 
measured  in  the  present  study  or  reported  by  others. 

©  2011  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  have  attracted 
much  attention  in  the  past  decade  as  a  promising  candidate  for 
clean  power  sources  in  automotive,  electronic,  portables,  and  sta¬ 
tionary  applications  [1].  Challenges  facing  the  commercialization 
of  PEMFCs  are  reliability,  cost,  power  density,  performance,  and 
membrane  lifetime.  Membrane  electrode  assembly  (MEA),  which 
is  the  heart  of  a  PEMFC,  is  composed  of  a  membrane  covered  by  cat¬ 
alyst  particles  sandwiched  between  two  porous  gas  diffusion  layers 
(GDLs)  on  the  sides  [1]. 

During  operation,  hydrogen  and  oxygen  are  consumed  in  the 
anode  and  cathode  sides  of  the  MEA,  respectively;  and  heat  and 
water  are  produced  in  the  cathode  side.  As  a  result  of  the  electro¬ 
chemical  reactions,  an  electrical  current  is  generated  from  anode  to 
cathode  via  an  external  circuit  that  can  be  used  for  various  appli¬ 
cations.  To  ensure  the  optimum  operating  condition  of  a  PEMFC, 
reactants  should  continuously  be  transported  towards  the  cata¬ 
lyst  layer,  where  the  electrochemical  reactions  occur.  The  excess 
produced  water  and  heat  should  be  removed  effectively  and  elec¬ 
trons  should  pass  from  the  bipolar  plate  towards  the  catalyst  layer 
or  in  the  reverse  direction  in  cathode  and  anode  sides,  respec¬ 
tively  [2].  The  GDL  plays  a  crucial  role  in  all  the  aforementioned 
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transports;  moreover,  it  provides  mechanical  support  for  the  MEA. 
As  a  result,  the  transport  properties  of  GDLs  are  extremely  impor¬ 
tant  in  assessing  the  overall  performance  of  PEMFCs.  The  mass 
transport  loss,  related  to  GDL  permeability,  is  the  major  cause  for 
limiting  the  maximum  power  density  of  a  PEM  fuel  cell.  As  such, 
an  in-depth  knowledge  of  permeability  of  GDLs  is  important  for 
optimizing  PEMFCs  operating  condition  [3-10]. 

In  a  PEM  fuel  cell  stack,  the  cell  components  are  compressed 
together  to  prevent  gas  leakage  from  the  system.  The  applied  pres¬ 
sure  changes  the  GDL  thickness  and  microstructure;  this  affects  the 
transport  properties  of  the  GDL  to  a  considerable  extent.  The  gas 
permeability  reduces  as  a  result  of  the  GDL  compression  while  the 
thermal  and  electrical  conductivities  increase  [5].  Therefore,  a  trade 
off  exists  in  the  design  and  optimization  process  of  MEAs  [9,11]. 

Water  management  is  one  of  the  current  challenges  in  increas¬ 
ing  power  density  of  PEMFCs.  Use  of  hydrophilic  GDLs  leads  to 
membrane  dehydration  and  high  saturation  which  causes  the  pores 
blockage  and  hindering  of  the  reactant  access  to  catalyst  layers 
[12-14].  To  prevent  these  issues,  the  GDL  is  commonly  treated 
with  hydrophobic  Polytetrafluoroethylene  (PTFE);  this  affects  the 
porosity  and  microstructure,  which  consequently  changes  the  gas 
permeability  of  the  treated  GDL  and  the  mass  transfer  rate  [15]. 
Thus,  model(s)  that  can  predict  the  permeability  of  GDLs  as  a 
function  of  porosity,  mechanical  compression,  and  PTFE  content 
is  important  for  designing  optimized  PEMFCs. 

Generally,  GDLs  are  porous  media  made  of  carbon  fibers  with 
diameter  of  7-10  p>m  in  the  form  of  carbon  papers  or  carbon  cloths 
[16].  These  materials  have  a  complex  microstructure  with  a  random 
distribution  of  fibers  which  makes  developing  analytical  solutions 
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Nomenclature 

d  Fiber  diameter,  m 

E(.)  Uncertainty  associated  with  the  measurement  of  a 
parameter 

GDL  Gas  diffusion  layer 

I<  Permeability,  m2 

P  Pressure,  Pa 

PEMFC  Polymer  electrolyte  membrane  fuel  cell 
PTFE  Polytetrafluoroethylene 

Q.  Volumetric  flow  rate,  m3  s-1 

S  Distance  between  centers  of  adjacent  fibers,  m 

t  Thickness  of  the  samples,  m 

UD  Volume-averaged  superficial  velocity,  m  s-1 

Greek  symbols 

<5min  The  minimum  distance  between  adjacent  fibers  sur¬ 
faces,  m 
s  Porosity 

li  Fluid  viscosity,  Ns  m-2 

ip  Solid  volume  fraction,  cp  =  1  -  s 

x  Tortuosity  factor 

co  PTFE  content 

Subscripts 

comp  Compressed  GDL 

0  Uncompressed,  non-treated  GDL 


for  flow-field  and  permeability  highly  unlikely.  Therefore,  the  per¬ 
meability  of  GDLs  has  been  mostly  studied  either  experimentally 
or  numerically  in  the  literature. 

There  are  few  experimental  studies  available  in  the  open  lit¬ 
erature  in  which  the  permeability  of  GDLs  has  been  investigated 
[2,15-20].  Ihonen  et  al.  [15]  showed  that  a  reverse  relation¬ 
ship  exists  between  the  in-plane  permeability  and  compression. 
Recently,  Leser  et  al.  [2]  studied  the  effects  of  mechanical  compres¬ 
sion  and  reported  the  in-plane  gas  permeability  as  a  function  of 
porosity  for  a  carbon  cloth,  a  non-woven  carbon  fiber  GDL,  and  a 
carbon  paper.  In  a  similar  work,  Gostick  et  al.  [  1 6]  measured  the  per¬ 
meability  of  several  commercial  GDLs  under  various  compressive 
loads  and  reported  the  in-plane  permeability  as  a  function  of  poros¬ 
ity.  They  also  reported  through-plane  permeability  of  several  GDLs. 
Recently,  Tahler  et  al.  [20]  and  Becker  et  al.  [19]  have  measured  the 
through  plane  permeability  of  compressed  Toray  carbon  papers. 
Gurau  et  al.  [18]  measured  in-plane  and  through-plane  permeabil¬ 
ities  of  a  carbon  paper  with  various  PTLE  contents  and  reported  a 
direct  relationship  between  PTLE  content  and  permeability.  On  the 
contrary,  the  experimental  observation  of  Park  et  al.  [21  ]  showed  a 
reduction  of  permeability  by  increasing  PTLE  content. 

On  the  theoretical  side,  Lattice  Boltzmann  simulations  of  gas 
flow  through  several  random  fibrous  structures  were  carried  out 
by  Vandoormaal  and  Pharoah  [22]  over  the  porosity  range  of 
0.6  <  s  <  0.8.  Their  numerical  results  suggested  a  significant  impact 
of  fiber  orientations  on  the  permeability  even  in  a  constant  poros¬ 
ity.  Based  on  the  numerical  results  for  structures  with  mixed 
fiber  orientations,  they  proposed  two  correlations  for  in-plane  and 
through-plane  permeabilities  of  the  considered  structures;  how¬ 
ever,  these  correlations  are  not  applicable  to  compressed  or  PTLE 
treated  GDLs.  Becker  et  al.  [  1 9]  used  3D  tomography  to  reconstruct  a 
GDL  and  a  numerically  efficient  pore  morphology  method  to  deter¬ 
mine  phase  distributions  and  to  deduce  permeability,  diffusivity 
and  thermal  conductivity  as  a  function  of  the  saturation  under  dif¬ 
ferent  compressive  loads.  Hao  and  Cheng  [23 ]  reconstructed  carbon 
papers  with  different  PTLE  contents  using  a  stochastic  method  and 


solved  the  pore-level  flow  by  Lattice  Boltzmann  method.  Shi  et  al. 
[24]  employed  the  fractal  theory  and  developed  models  for  per¬ 
meability  that  accounted  for  the  microstructures  of  GDLs  in  terms 
of  two  fractal  dimensions  and  proposed  a  relationship  for  perme¬ 
ability  as  a  function  of  tortuosity,  fractal  dimensions,  pore  area, 
pore  size  distribution,  and  effective  porosity.  Their  model,  however, 
requires  several  geometrical  parameters  (e.g.  fractal  parameters) 
that  should  be  known  beforehand.  A  similar  approach  was  also 
adopted  by  He  et  al.  [25].  Recently,  Tamayol  and  Bahrami  [26]  have 
proposed  a  blending  technique  for  predicting  the  in-plane  perme¬ 
ability  of  GDLs.  They  treated  the  medium  as  a  mixture  of  fibers 
parallel  and  perpendicular  to  the  flow  direction  and  estimated  the 
permeability  as  a  blend  of  the  parallel  and  normal  permeability  of 
unidirectionally  aligned  fibers  [27,28]. 

Through-plane  permeability  of  uncompressed  carbon  papers 
can  be  estimated  using  the  existing  models  for  the  permeability  of 
planar  fibrous  structures.  Based  on  the  analogy  between  electrical 
and  flow  conductions,  Tomadakis  and  Robertson  [29]  developed 
a  model  for  permeability  of  randomly  distributed  overlapping 
fibers  in  composite  reinforcements.  Recently,  Tamayol  and  Bahrami 
[30]  used  a  scale  analysis  technique  and  related  the  permeabil¬ 
ity  of  fibrous  media  to  the  microstructure  geometrical  parameters 
including  tortuosity,  pore  diameter,  and  porosity.  They  have  suc¬ 
cessfully  compared  their  model  with  experimental  data  collected 
from  various  sources  for  a  variety  of  materials.  However,  the  models 
of  [29,30]  do  not  include  the  effect  of  PTLE  content  and  mechanical 
compression  on  the  through-plane  permeability. 

Therefore,  there  is  a  need  for  a  general  model  that  can  accu¬ 
rately  predict  the  through-plane  permeability  of  GDLs  as  a  function 
of  compression  ratio  and  PTLE  content.  In  addition,  more  experi¬ 
mental  investigations  are  required  to  verify  the  developed  model 
and  determine  effects  of  mechanical  compression  and  PTLE  con¬ 
tent  on  the  through-plane  permeability  of  GDLS.  In  this  study,  the 
through-plane  permeability  of  several  compressed  GDLs  is  mea¬ 
sured  experimentally.  In  addition,  effect  of  PTLE  content  on  the 
permeability  is  investigated.  The  model  of  Tamayol  and  Bahrami 
[30]  is  modified  to  accommodate  the  effects  of  mechanical  com¬ 
pression  and  PTLE  contents. 


2.  Experimental  approach 

2.1.  Tested  samples 

The  tested  GDLs  were  purchased  as  Teflon  treated  carbon  papers 
from  different  manufacturers.  As  a  result  of  their  similar  produc¬ 
tion  procedure,  despite  the  different  thicknesses,  TGP-H-90  and 
120  have  similar  microstructures.  TGP-H-120  was  obtained  with 
0,  5,  20,  and  30%  PTLE  content  while  the  tested  TGP-H-90  have  a 
0%  PTLE  content.  Lig.  la  shows  a  SEM  image  of  TGP-H-120  with  5% 
PTLE  content.  The  SGL  Sigracet  1 0AA,  shown  in  Lig.  1  b,  was  obtained 
with  0%  PTLE  content. 

The  porosity  of  the  samples  which  was  reported  by  manufac¬ 
turer  and  Gostick  et  al.  [16]  are  used  in  the  present  study.  Liber 
diameter  was  reported  by  Gostick  et  al.  [16]  and  was  also  inde¬ 
pendently  calculated  using  scanning  electron  microscopy  (SEM) 
images.  The  thickness  of  the  samples  was  measured  using  a 
micrometer;  the  maximum  uncertainty  associated  with  the  mea¬ 
surements  is  less  than  lOp^m.  The  properties  of  the  purchased 
samples  are  listed  in  Table  1. 

2.2.  Test  apparatus 

An  air  permeability  test  bed  was  designed  and  built  for  measur¬ 
ing  the  through-plane  permeability,  see  Lig.  2.  The  high  pressure 
air  was  supplied  by  a  high  pressure  air  tank  and  controlled  using  a 
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Table  1 

Physical  properties  of  different  carbon  papers  used  in  the  present  study:  porosity  and  PTFE  content  provided  by  manufacturer;  the  measured  thicknesses  and  fiber  diameters. 


GDL 

Thickness  (pirn) 

Porosity3 

Fiber  diameter  (p,m) 

PTFE  content  (%) 

SGL-Sigracet  10AA 

390-410 

0.88 

9.2 

0 

TGP-H-90 

255-270 

0.79 

9.3 

0-30 

TGP-H-120 

370-385 

0.79 

9.3 

0 

a  The  porosity  is  reported  for  GDLs  with  0%  PTFE  content. 


digital  air  pressure  regulator;  the  output  pressure  of  the  air  regu¬ 
lator  valve  was  set  to  6  psi  during  all  experiments.  Two  aluminum 
blocks  were  machined  and  drilled  with  a  25  mm  diameter  hole. 
The  blocks  were  used  as  sample  holders  and  GDL  samples  were 


sandwiched  between  the  blocks.  An  O-ring  was  used  between  the 
two  aluminum  blocks  to  prevent  air  leakage  from  the  test  section. 
The  up  and  down  stream  air  flow  were  connected  to  a  differen¬ 
tial  pressure  transducer,  PX277-01 D5V,  supplied  by  Omega  (Omega 
Inc.,  Laval,  QC,  Canada).  The  pressure  transducer  was  connected  to 
a  PC,  using  a  DAQ  purchased  from  Omega  (Omega  Inc.,  Laval,  QC, 
Canada),  where  the  pressure  drop  values  were  recorded.  The  accu¬ 
racy  of  the  pressure  transducer  was  1%  of  its  full  scale  measuring 
range  (0-1  inch  of  water).  The  air  flow  rate  passing  through  tested 
GDLs  was  measured  using  a  flow  meter  provided  by  Omega  (Omega 
Inc.,  Laval,  QC,  Canada).  The  accuracy  and  the  range  of  the  air  flow 
meter  were  3%  and  0-10 1pm,  respectively. 

The  carbon  papers  were  cut  into  circular  samples  of  3  cm  diame¬ 
ter  for  the  experiments.  To  investigate  the  compression  effects,  GDL 
samples  were  compressed  using  a  GUNT  WP  300  Universal  Material 
Tester  for  1 5  min;  the  maximum  applied  compressive  force  was  less 
than  20 1<N  (30  MPa).  The  thickness  of  the  compressed  samples  was 
measured,  using  the  same  micrometer,  after  the  load  was  released 
to  determine  the  variation  in  its  thickness. 


2.3.  Data  analysis 


Gas  permeability,  I<,  is  defined  using  Darcy  equation  [31  ]: 


A  P 

~T 


(1) 


where  UD  is  the  volume  averaged  velocity  through  porous  media, 
fi  is  viscosity,  A P  is  the  pressure  drop  across  the  sample,  and  t 
is  the  GDLs  thickness.  Eq.  (1)  is  valid  for  incompressible,  steady, 


Fig.  1.  SEM  images:  (a)  TGP  120  with  5%  PTFE  content;  (b)  SGL  Sigracet  10AA  and  Fig.  2.  The  air  permeability  testbed:  (a)  schematic  of  the  apparatus  (exploded  view) 

(c)  compressed  TGP  1 20  with  5%  PTFE  content.  and  (b)  actual  test  setup. 
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Fig.  3.  Measured  pressure  drops  for  samples  of  compressed  TGP-H-120. 


constant  properties,  single-phase  (no-surface  tension  forces),  and 
low  Reynolds  number  flows.  The  Reynolds  number  based  on  the 
fiber  diameter,  Re  =  pUDd//i ,  was  lower  than  0.1  in  the  present 
study;  therefore,  the  Darcy  equation  was  applicable  for  analyzing 
the  experimental  results.  The  maximum  measured  pressure  drop 
across  the  samples  was  less  than  200  Pa.  As  such,  it  was  assumed 
that  the  compressibility  effects  for  such  small  values  of  pressure 
drop  can  be  negligible  and  the  Darcy  equation  in  the  form  of  Eq.  ( 1 ) 
was  applicable  in  the  present  study.  UD  was  calculated  by  dividing 
the  volumetric  flow  rate  by  the  sample  cross-sectional  area. 

Assuming  Darcy’s  law  in  a  porous  structure  implies  a  linear  rela¬ 
tionship  between  the  pressure  drop  and  the  fluid  velocity  in  the 
media.  As  an  example,  this  linear  relationship  can  be  observed  in 
the  experimental  results  for  samples  of  compressed  TGP-H-120  in 
Fig.  3  which  justifies  Darcy’s  assumptions. 

The  measured  values  of  pressure  drop,  A Ptota/,  during  the  exper¬ 
iment  were: 


AP total 


A  P, 


sample 


+  AP, 


min 


(2) 


where  A Psampie  is  the  pressure  drop  associated  within  the  sam¬ 
ple  and  APmin  is  the  pressure  drop  related  to  minor  losses  in  the 
testbed,  i.e.,  exit  and  entrance  section  of  the  sample  holder.  To 
account  for  the  minor  pressure  losses,  the  pressure  drop  in  the  test 
section  was  measured  once  without  any  GDL  samples  to  determine 
APmin  for  the  entire  range  of  flow  rate  considered  in  the  present 
experimental  work. 

The  uncertainty  associated  with  the  through-plane  permeabil¬ 
ity,  calculated  based  on  the  measured  variables  using  the  Darcy’s 
equation,  can  be  estimated  from  [32]: 


E(I<) 

I< 


+ 


(E(A  P) 
V  AP 


(3) 


Table  2 

Uncertainty  values  for  measured  parameters. 


Parameter 

Uncertainty 

t 

10  |xm 

AP 

1%  of  full  scale 

Q 

3%  of  full  scale 

£ 

5% 

I< 

9% 

3.  Theoretical  model 


Tamayol  and  Bahrami  (TB)  [30]  employed  a  scale  analysis  and 
offered  a  scale  for  permeability  of  planar  fibrous  microstructure, 
shown  in  Fig.  4,  as: 


KsS2  •  t 

6°min L 


(4) 


where  e  was  the  porosity,  r  was  the  tortuosity  factor,  and  <5min  was 
half  of  the  minimum  opening  between  two  adjacent  fibers  (S  -  d), 
see  Fig.  4.  They  estimated  the  tortuosity  factor  from  the  Bruggeman 
equation  [33]: 


(5) 


where  <p  was  the  solid  volume  fraction,  i.e.,  <p  =  1  -  s.  Substituting  for 
geometrical  parameters,  the  following  relationship  was  reported 
for  the  dimensionless  through-plane  permeability: 


T.  =  0.008^(1  —  f) 


(6) 


The  constant  value  in  Eq.  (6),  i.e.,  0.008,  was  found  through  com¬ 
parison  with  experimental  data  collected  from  different  sources  for 
permeability  of  various  planar  fibrous  materials.  Hao  and  Cheng 
[23]  numerically  calculated  the  tortuosity  factor  for  carbon  papers 
and  proposed  the  following  correlation: 


r  =  1  +0.72 


1 


(£  —  0.11) 


0.54 


(7) 


Our  analysis  showed  that  if  the  constant  value  in  Eq.  (6)  is  replaced 
by  0.012  and  r  is  calculated  from  Eq.  (7)  the  resulting  equation: 


T.  =  0.012(1  —  cp) 


(  7T  V  n  71  -1 

\4  (fi)  4<p  + 


1+0.72 


(0.89  —  (p) 


0.54 


(8) 


can  predict  the  experimental  data  for  GDLs  more  accurately. 

To  enable  the  model  of  Tamayol  and  Bahrami  [30]  to  accom¬ 
modate  the  effects  of  compression  factor  and  PTFE  content,  the 


where  E(.)  is  the  uncertainty  in  the  measurement  of  each  parame¬ 
ter;  these  values  are  listed  in  Table  2  The  maximum  uncertainty  in 
the  calculated  values  of  permeability  is  estimated  to  be  9%. 


Fig.  4.  The  2D  unit  cell  considered  in  the  present  study. 
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Fig.  5.  Effect  of  compression  ratio,  tcom(J/t0  on  the  variation  of  permeability.  FiS- 7-  Effect  of  PTEE  content  on  the  through-plane  permeability  of  two  set  of  TGP- 

H-l  20  samples  with  various  PTFE  contents. 


relationship  between  these  properties  and  the  solid  volume  frac¬ 
tion  should  be  determined. 

Fig.  lc  shows  that  mechanical  compression  does  not  change 
the  shape  of  fibers.  As  such,  it  can  be  assumed  that  during  the 
compression  process  only  the  thickness  of  the  original  fibrous  sam¬ 
ples  changes  while  the  volume  of  the  solid  carbon  fiber  remains 
constant.  Therefore,  the  relationship  between  the  solid  volume 
fraction,  (pCOmp ,  of  a  compressed  sample  with  the  original  value  for 
an  uncompressed  carbon  paper,  (p0,  can  be  expressed  as  [23]: 


[ Pcomp  —  <Po 


Jo_ 

tcomp 


(9) 


filled.  The  final  porosity,  Sptfe,  can  be  expressed  approximately  as 
a  function  of  PTFE  content  oo  [23]: 


SpTFE  =  So  —  a 


co(  1  -£p) 

1  -  CO 


(10) 


where  so  is  the  original  porosity  before  PTFE  treatment.  Hao  and 
Cheng  [23]  suggested  a  =  0.9  as  the  density  ratio  of  the  carbon 
fiber  and  the  PTFE  [34].  Employing  Eqs.  (9)  and  (10),  one  can  pre¬ 
dict  effects  of  PTFE  and  compression  on  the  permeability  on  the 
through-plane  permeability  of  GDLs  from  TB’s  model. 


4.  Results  and  discussions 


where  tcomp  and  t0  are  the  compressed  and  uncompressed  GDL 
thicknesses,  respectively.  If  the  PTFE  is  added  on  the  carbon  paper 
GDL,  the  pore  volume  is  randomly  filled  by  the  PTFE.  It  is  postulated 
that  PTFE  changes  the  porosity  and  some  pores  in  the  medium  are 


The  permeability  of  various  tested  samples  is  calculated  using 
Darcy’s  law,  Eq.  (1).  The  porosity  of  the  compressed  samples 
and  GDLs  treated  with  PTFE  are  calculated  using  Eqs.  (9)  and 
(10),  respectively.  The  through-plane  permeability  of  uncom¬ 
pressed  TGP-H-90  was  measured  as  8.7  x  10~8  m2  which  is  in  good 


Fig.  6.  Comparison  of  the  proposed  model  with  the  experimental  data  for  com¬ 
pressed  GDLs  measured  in  the  present  study  or  collected  from  various  sources. 


Fig.  8.  Comparison  of  the  proposed  model  with  the  experimental  data  for  TGP-H- 
120  with  various  PTFE  contents. 
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agreement  with  the  value  of  8.99  x  10-8  m2  reported  by  Gostick 
et  al.  [16]. 

The  effect  of  mechanical  compression  and  variation  of  GDL 
thickness  on  the  permeability  is  shown  in  Fig.  5.  It  can  be  seen 
that  there  is  a  linear  relationships  between  the  ratio  of  compressed 
to  uncompressed  permeability  of  the  measured  GDLs  and  the 
compression  ratio,  tCOmp/fo-  The  experimental  data  for  the  through- 
plane  permeability  of  compressed  GDLs,  from  the  present  study  or 
reported  by  others,  are  plotted  in  Fig.  6  and  compared  with  the 
present  model  (TB  model),  Eq.  (8),  and  the  correlation  proposed  by 
VanDoormaal  and  Pharoah  [22]  for  the  through-plane  permeabil¬ 
ity  of  GDLs.  It  can  be  seen  that  the  TB  model,  Eq.  (8),  has  a  better 
agreements  with  the  experimental  data  for  compressed  GDLs. 

The  effect  of  PTFE  content  on  the  through-plane  permeability 
of  GDLs  is  presented  in  Fig.  7.  It  can  be  seen  that  a  reverse  rela¬ 
tionship  exists  between  the  PTFE  content  and  the  through-plane 
permeability.  The  experimental  data  for  TGP-H-120  with  various 
PTFE  contents  are  plotted  in  Fig.  8  and  compared  with  TB  model. 
Porosities  of  the  samples  are  estimated  from  Eq.  (10).  The  perme¬ 
ability  of  the  PTFE  treated  GDLs  has  a  reverse  relationship  with  the 
PTFE  content.  It  can  be  seen  that  Eq.  (8)  captures  the  trends  of  the 
data.  Overall,  it  can  be  concluded  that  the  modified  TB  model  can 
be  used  in  design  and  optimization  of  PEMFCs. 

5.  Conclusions 

An  experimental  study  has  been  carried  out  to  investigate 
the  effects  of  compression  and  PTFE  content  on  the  single-phase 
through-plane  permeability  of  GDLs  of  PEMFCs.  A  custom-made 
testbed  was  built  and  several  GDL  samples  were  obtained  from  var¬ 
ious  manufacturers  and  tested  over  a  wide  range  of  compression 
ratios.  Moreover,  the  through-plane  permeability  of  TGP-H-120 
carbon  papers  with  various  PTFE  contents  was  measured.  A  scale 
analysis  model  originally  proposed  by  Tamayol  and  Bahrami  [30] 
was  modified  for  predicting  the  permeability  of  the  PEMFC  GDLs. 
The  present  model  takes  into  account  the  PTFE  content  and  com¬ 
pression  effects  of  the  permeability. 

The  highlights  of  the  present  study  can  be  summarized  as: 

•  A  reverse  relationship  exists  between  the  through-plane  perme¬ 
ability  of  GDLs  and  both  the  compression  ratio  and  PTFE  content. 

•  The  revised  model  of  Tamayol  and  Bahrami,  Eq.  (8)  captures  the 
trends  observed  in  experimental  data  for  permeability  of  com¬ 
pressed  GDLs. 

The  reported  experimental  data  and  the  proposed  model  can  be 
used  to  guide  the  design  and  optimization  of  PEMFCs,  and  can  be 
readily  implemented  into  fuel  cell  models  that  require  specification 
of  the  through-plane  gas  permeability  of  the  GDL. 
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